Semi-classical magnetoresistance in weakly modulated magnetic fields 



O 

o 
o 

a 

in 

(N 



X3 

o 



> 

in 
m 

o 
o 
o 



I 

o 
o 



X 



A. MatuliJ and F. M. Peeterfi 
Department of Physics, University of Antwerp (UIA), B-2610 Antwerp, Belgium, 

(February 1, 2008) 



The semi-classical conductance of a two-dimensional elec- 
tron gas is calculated in the presence of a one-dimensional 
modulated magnetic field with zero average. In the limit of 
small magnetic field amplitudes (B) the contribution of the 
magnetic modulation to the magnetoresistance increases as 
B^^^ in the diffusive limit, while the increase is linear in B 
in the ballistic regime. Temperature does not influence the 
power law behavior but it decreases the prefactor of this func- 
tional behavior. 

PACS numbers: 73.23.-b, 73.40. Gk, 75.60.-d 



In recent years there has been an increased interest |^] 
in hybrid systems which promise to increase the function- 
ality of present day semiconductor devices. One example 
of such type of systems are those in which semiconductors 
and magnetic materials are combined where the magnetic 
material provides a local magnetic field which influences 
locally the motion of the electrons in the semiconductor. 
The latter is usually a heterostructure which contains a 
two-dimensional electron gas (2DEG). The 2DEG acts as 
a detector measuring the magnetic state of the magnetic 
material. Previously, the coupling between such a non 
homogeneous magnetic field and the 2DEG was demon- 
strated lH in the case the magnetic field {B) was directed 
perpendicular to the 2D plane. In this case one has a 
modulation of the i?-field on top of a homogeneous back- 
ground field and the influence of the _B-field modulation 
on the 2DEG is relatively weak . 

When the magnetic field is directed parallel to the 
2DEG the magnetic material is magnetized parallel to 
the 2D plane which leads to fringing fields near the edge 
of the magnetic material having a non zero magnetic field 
component perpendicular to the 2D layer. Those fields 
form a magnetic barrier for the electron motion in 
the 2D plane. Because now there is no background per- 
pendicular magnetic field for the 2DEG the influence on 
the resistance of such magnetic barriers is much more 
pronounced and large increases in the magnetoresistance 
have been found 

In the present work we investigate the magnetotrans- 
port in weak modulated magnetic fields for which the 
average B-field is zero. For the case the typical magnetic 
energy is much smaller than the Fermi energy a semi- 
classical analysis is applicable. We find that in the diffu- 
sive regime the correction to the magnetoresistance ex- 
hibits a non analytical behavior in the limit of small mag- 
netic field amplitudes which differs from the behaviour in 



the ballistic regime. 

We consider electrons moving in a two-dimensional 
(2D) xy-plane. The magnetic field, directed along 
the z-direction, is periodic along the si-direction B — 
5(0, 0, 5(ai/Zo)) with period where h{x) is a periodic 
{b{x -|- 1) = h{x)) dimensionless function describing the 
magnetic field modulation with zero average value. 

In a semi-classical analysis the electron motion in a 
magnetic field is described by the following Hamiltonian 
(or its energy): 
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where m is the electron effective mass, p — {px,Py) is 
the electron canonical momentum, and the dimensionless 
periodic function 
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characterizes the vector potential A = Blo{0, a{x/lo),0) 
which is taken in the Landau gauge. The quantum energy 
spectrum of electrons in modulated magnetic fields was 
studied in Refs. §,|l§. 

We restrict our analysis to the case where the electron 
transition through a single period is ballistic, i. e. the 
mean free path I = vpT ^ Iq {vp is the electron Fermi 
velocity, and t the relaxation time) and the motion in 
the sample is diffusive, i.e. the mean free path is smaller 
than the size of the sample, [l <C Lx,Ly). For diffusive 
transport and in the limit of small magnetic fields {oJcT <C 
1) the expression for the average conductivity tensor is 
given by the following integral over the electron phase 
space {x,px,py) 
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where the symbol fpi^) — {exp{(/i — £)/kT} + 1}^^ 
stands for the equilibrium electron Fermi-Dirac distribu- 
tion function with T the temperature and /x the chemical 
potential which equals the Fermi energy Ep = Trfi^n/m 
in the zero temperature limit, where n is the 2D electron 
density. Note that the coordinate y is excluded from the 
phase space as the system is homogeneous in that direc- 
tion. The expressions for the electron velocities follow 
from the Hamilton equations of motion 
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First, let us consider the conductivity in the zero tem- 
perature Umit when the derivative of the Fermi function 
reduces to a ^-function. The component ayy can be cal- 
culated straightforwardly, as the sample is homogeneuos 
along the y-direction all trajectories have to be taken into 
account. Inserting expression into the conductivity 
expression we obtain 



where the symbol a^°' stands for the snake orbit con- 
tribution. The latter term defines the decrement of the 
conductivity due to the modulated magnetic field, and 
in the limiting case of small magnetic fields it is propor- 
tional to the increase of the magnetoresistance due to the 
magnetic field modulation. 

The snake orbits are located above the plane A and 
below the plane B. Those planes are defined by Py{A) — 
^2mEF + [eBla / c)amin, and Py{B) = — \/2to£'f -I- 
(ei?Zo/c)amax, with Omin = min{a(x)} and amax = 
max{a(x)} the extremal points of the vector potential 
a{x). Thus we can write 
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which leads to the well-known magnetic field independent 
result 



(6) 



We assumed that the relaxation time depends only on 
the electron energy and therefore for the case of zero 
temperature it can be replaced by the constant value r = 
t{Ep). Thus the weak magnetic field modulation in the 
x-direction does not change the conductivity in the y- 
direction, as expected. 

The calculation of axx is more complicated because for 
sufhciently strong magnetic fields (or small electron ve- 
locity, but such that TluJc <C Ep) some electrons can be 
forced into snake orbits. The electron motion on such or- 
bits oscillates in the x-direction around the average value 
xq. Therefore, such electrons do not contribute to jx, and 
consequently, in the expression for the conductivity Gxx 
those snake orbits have to be excluded. The classification 
of all possible electron orbits are given in Fig. |l| where the 
Fermi surface e{px,Py,x) = Ep is plotted. In the case of 
zero temperature only electrons with trajectories on the 
Fermi surface contribute to the conductivity integral. As 
the energy (|l|) does not depend on y, the momentum pj, is 
conserved, and consequently the trajectories are defined 
by the intersection of the above Fermi surface with the 
Py = C*'^ planes. It is apparent from Fig. |^ that there are 
two types of trajectories. The trajectories as indicated 
by symbol D are able to run along the whole a;-axis. 
Such electrons are moving along open trajectories and 
they will contribute to the current along the x-direction. 
The other trajectories, indicated by symbol E, are closed 
and they correspond to the snake orbits. Thus we can 
separate the conductivity into two parts 
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where the contribution from the snake orbits above the 
A plane is given by 
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Here the angular integration has to be performed over 
the white sections of the circle C in Fig. or equiva- 
lently, along the bold arch shown in Fig. pT The limit- 
ing angle is defined as (po = arccos(I — Apy/y/2mEp), 
with Apy — (eBlo/c) {a(x/lo) — flmin}- In the asymp- 
totic case of small magnetic field modulations it becomes 



ipo ^ y 2Apy/ \J2raEp. In the latter case the integration 
in expression (0) leads to 
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where Bo = {2t: c^/rnEp) / {elo), and 
8-21/4^ 
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The integration over the snake orbits below the plane B 
leads to the same expression (|l^) except that now the 
coefficient ca has to be replaced by 
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Thus the resistance change due to the modulated mag- 
netic field becomes ARxx/Ro = {ca + cb){B / BoY^^ 
where i?o is the resistance in the absence of any mag- 
netic field modulation. 
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As a special case let us consider a simple cosine peri- 
odic magnetic field modulation (with period Iq) b{x) = 
cos(27ra;), which leads to a{x) — sin(27rx)/27r. In this 
case the coefficients can be easily evaluated 
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In the present classical ballistic situation an electron 
which passed through the first magnetic barrier will also 
pass through the other barriers. As a consequence the 
above result can also be applied to the one barrier sit- 
uation. The periodic oscillating Fermi surface shown in 
Fig. d reduces now to a single step. The trajectories of D- 
type (i.e. the open orbits) contribute to the a^x conduc- 
tivity, but the snake orbits (trajectories of A-type) are 
replaced by trajectories which reflect from the barrier. 
The integral in expression ( p^ ) must now be evaluated 
over those reflected trajectories. When the barrier thick- 
ness lo <C Lx, the integral j^'°J'^i^dx{a{x /lo) — flmin)'^^^ 
becomes the sample length multiplied by the total vector 
potential increment over the barrier: irciflmax — amin}'^^^, 
where Omax = a{Lx/2) and a,nin = a(-Lx/2). 

A single magnetic field barrier in a 2DEG, created ex- 
perimentally by parallel magnetization of magnetic strips 
placed on top of the 2DEG, can be represented by [p]Jll| 
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The magnetic strip has a thickness D with magnetization 
M and is placed a distance d from the 2D electron sys- 
tem. The vector potential is obtained by integrating (14) 
over {~Lx/2, Lx/2) which gives Amax-^min = fJ.oMD/2 
where use has been made of ^ d, D which is valid in 
typical experimental situations and A — Bloa. Finally, 
we obtain for the contribution of the reflected trajectories 
to the conductivity 
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The main feature of the obtained magnitoresistance 
is its non analytical behavior 5^/2 fQp small magnetic 
field amplitudes. It is remarkable that it does not de- 
pend on the actual form of the modulating field, but it 
is determined by the density of snake orbits (or reflected 
trajectories) at the Fermi surface. 

For non zero temperature, the previous effect will be 
suppressed by thermal fluctuations. To generalize our 
results to non zero temperature we have to replace any 
function G{Ep) which depends on the Fermi energy Ep 
by the corresponding average over the derivative of the 
Fermi function 
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Consequently, taking into account cTyy (||) and the depen- 
dence of Bo on Ep we obtain for the snake orbit contri- 
bution to the conductivity 
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where we assumed that the relaxation time does not de- 
pend on the energy. In the limit of small temperatures 
the integral can be evaluated analytically and we arrive 
at the final expression for the conductivity along the di- 
rection of the magnetic field modulation 
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Notice that temperature decreases the nonanaliticity co- 
efficient Co (with about 30% when fc^T = /i) but does 
not influence the power law dependence. 

Tunneling through magnetic barriers in the ballistic 
regime was studied numerically in Refs. QJ^. For weak 
barrier only electrons impinging on the magnetic barrier 
under an angle 9 « tt/2 (i.e. cos(6') w (7''/2 — 9) ~ cf)) 
are reflected and thus we obtain from Eq. (6) of Rcf. § 
the conductance change due to tunneling through a weak 
magnetic barrier 
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where Go = e^mvpLy/h? . Notice that the difference 
with the diffusive case is the power of the angle in the 
expression for the conductance/conductivity. In the case 
of ballistic transport the conductance is proportional to 
Vx, while for diffusive transport, see Eq. (|^), we have 
1)2 . This difference leads to a linear i?-dependence in the 
ballistic regime 
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The above expression is obtained for the single barrier 
(|l|) and agrees with Eq. (3) of Ref. §. 

Lets compare these results with the experiments of 
Refs. [§§. In Ref. § a Co strip of thickness D = QOnm 
was placed a distance d = 35nm from a 2DEG formed 
in a GaAs-heterostructure with Ep — 15.7meV. The 
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magnetization of the Co strip ||T^ was fioM « 9B with 
a saturation magnetization of 1.6T. In this experiment 
one is in the diffusive regime and inserting these val- 
ues in Eq. (15) we obtain AR^^/Rq = 4.3B(T)3/2. 
The estimated zero field resistance was Rq = 2.3557 
which results into AR.^^{n) = 10.1B(T)3/2 and agrees 
with the experimental low magnetic field behavior 
AR^^in) = QB{Tf/\ 

In contrast, the experiments of Vancura et al |^] 
are closer to the ballistic regime (they have been per- 
formed on shorter samples, i.e. Ly = 34/im). They 
placed rectangular Co dots of thickness D — lOOnm 
above a 2DEG of width Ly ^ W = 20nm in a GaAs- 
heterostructure with Ep — 19.7meV. Inserting these 
values in Eq. (19) we obtain AR^x/Rq — 3.6i?(T) where 
i?o = i/Go = l.Om. This gives AR^^i^) = 3.96B(T) 
which is a factor 2.2 smaller than the experimental result 
ARxxi^) = S.75B{T). At the saturation field of 1.5T we 
find theoretically AR^x = 5.94S7 which compares with 
the experimental result ARxx — 3.551 and which is a fac- 
tor 1.7 smaller. Thus on the average we find a reasonable 
agreement between theory and experiment but it is clear 
that all details are not yet fully understood. The dis- 
crepancy may be due to the fact that the Co dots are 
not homogeneously magnetized. Nevertheless, the linear 
magnetic field dependence is nicely reproduced. 

In conclusion we obtained the change in the magneto- 
resistance due to the presence of magnetic barriers with 
zero average magnetic field. We found that for diffusive 
transport the magneto-resistance increases as c{T)B^^^ 
with the amplitude of the magnetic barrier (B) where the 
coefficient c(T) decreases with increasing temperature. 
This result is different from the ballistic regime where 
the increase is linear in B. 
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FIG. 1. Fermi surface (thick solid curve) in phase space. 
Planes A and B are surfaces for the constant of motion py 
delimiting the open and closed orbits. Curve E (D) is an 
example of a closed (open) orbit. 



FIG. 2. Contour for angular integration in expression (|9|). 
Line A corresponds to the Fermi surface intersection indicated 
by plane A in Fig. 1. 
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